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Using the Monte Carlo simulation, we study the percolation and jamming of oriented linear k- 
mers on a square lattice that contains defects. The point defects with a concentration, d, are placed 
randomly and uniformly on the substrate before deposition of the k-mevs. The general case of 
unequal probabilities for orientation of depositing of k-mevs along different directions of the lattice 
is analyzed. Two different relaxation models of deposition that preserve the predetermined order 
parameter s are used. In the relaxation random sequential adsorption (RRSA) model, the deposition 
of k-mers is distributed over different sites on the substrate. In the single cluster relaxation (RSC) 
model, the single cluster grows by the random accumulation of k-mers on the boundary of the cluster 
(Eden-like model). For both models, a suppression of growth of the inhnite (percolation) cluster 
at some critical concentration of defects dc is observed. In the zero defect lattices, the jamming 
concentration pj (RRSA model) and the density of single clusters ps (RSC model) decrease with 
increasing length k-mers and with a decrease in the order parameter. For the RRSA model, the 
value of dc decreases for short k-mers {k < 16) as the value of s increases. For k = 16 and 32, 
the value of dc is almost independent of s. Moreover, for short k-mers, the percolation threshold is 
almost insensitive to the defect concentration for all values of s. For the RSC model, the growth 
of clusters with ellipse-like shapes is observed for non-zero values of s. The density of the clusters 
Ps at the critical concentration of defects dc depends in a complex manner on the values of s and 
k. An interesting Ending for disordered systems (s = 0) is that the value of ps tends towards zero 
in the limits of the very long k-mers, k ^ oo and very small critical concentrations dc 0. In this 
case, the introduction of defects results in a suppression of k-mer stacking and in the formation of 
‘empty’ or loose clusters with very low density. On the other hand, denser clusters are formed for 
ordered systems with ps ~ 0.065 at s = 0.5 and ps ~ 0.38 at s = I.O. 
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I. INTRODUCTION 

Deposition of large particles such as colloids, polymers 
or nanotubes on substrates can be considered and stud¬ 
ied as random sequential adsorption (RSA) [l|. During 
RSA, objects randomly deposit on the substrate; this 
process is irreversible, and the newly adsorbing objects 
cannot overlap or pass through previously deposited ones. 

The substrate may be prepatterned (e.g., see 0 ) or in¬ 
clude some impurities (defects) (e.g., see 0 ). The ad¬ 
sorbed objects may be identical or present a mixture 
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of objects of different sizes and shapes (e.g., see Hi). 
Moreover, anisotropy can be introduced by postulating 
unequal probabilities for the deposition of elongated ob¬ 
jects along different axes (e.g., see |,i). The anisotropy 
of deposition can reflect the influence of external fields, 
flows or anisotropy of the substrate. The adsorption of 
the elongated particles in the presence of external fields 
produces anisotropic layers (e.g., see 05 HI)- 

Very often, during RSA, a substrate can be considered 
as a discrete space, e.g., a regular or irregular lattice, or 
as a tree. The simplest, but most commonly used in¬ 
stance of a discrete substrate is the square lattice. The 
objects selected to be adsorbed onto a square lattice are 
generally linear in shape. Such linear segments (also de¬ 
noted as needles, linear chains, stiff rods or k-mers) con¬ 
sist of k successive connected sites. 

If deposition of the objects can continue indefinitely, a 
jamming state is reached. At the jamming state, there 
are still voids on the substrate between the previously de- 
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posited objects, but their sizes or shapes are not sufficient 
to allow the deposition of even one additional object. 

For a perfect lattice (a lattice without any defects), 
the jamming coverage, pj, is the fraction of the sites oc¬ 
cupied by deposited objects. For a diluted (disordered, 
disturbed) lattice (a lattice containing defects or impuri¬ 
ties), there are several ways to define the jamming cover¬ 
age [ol-lll] . In the present research, we have used only the 
pure object jamming limit and denoted it as the jamming 
concentration, p^, for short. The pure object jamming 
limit, pj, is defined as the fraction of the total number of 
lattice sites occupied by the deposited objects 0- 

Ben-Naim and Krapivsky obtained an important an¬ 
alytical result for the jamming concentration in a one¬ 
dimensional case 

If the concentration of the deposited objects on the 
substrate is sufficiently large, there may be a path 
through the objects from one side of the system to its 
opposite side. Below this concentration, there is no span¬ 
ning path through the system, while above this concen¬ 
tration, there is a connected component of the order of 
the size of the system. This concentration is called the 
percolation threshold. (See, e.g., [l2|for the details; art 
of the state may be found, e.g., in [i3-) For some sys¬ 
tems, percolation never occurs, even at the jamming con¬ 
centration (e.g., when the adsorbed layer is produced by 
the deposition of equally-sized squares on a square lattice 
and where the side of the square is greater than the length 
of 3 sites 0 ). In fact, the percolation threshold corre¬ 
sponds to a phase transition, e.g., insulator-conductor. 
Different definitions of the percolation threshold are used 
for lattices with defects similar to the definitions of the 
jamming concentration. In the present research, we have 
used the term ‘percolation threshold’ (pc) as the ratio of 
the sites occupied by the objects to the total number of 
lattice sites. 

Cornette et al. 00 investigated numerically the 
percolation of polyatomic species in the presence of im¬ 
purities on a square lattice with periodic boundary con¬ 
ditions. Bond and site percolation problems have been 
taken in consideration. Linear k-meis as well as so-called 
SAW /c-mers, i.e. segments of self avoiding walk, have 
been studied up to /c = 9. A phase diagram where the 
critical concentration of impurities is plotted as a func¬ 
tion of k has been offered. 

The kinetics of the random sequential adsorption of 
line k-meis (with values of k up to 64) has been studied 
on a disordered substrate occupied by the point impuri¬ 
ties 0 . The coverage of the surface and the jamming 
limits are calculated by the Monte Carlo method. The 
coverage has an asymptotically exponential behavior at 
low concentration of the impurities. The jamming limits 
depend on the concentrations of the impurities, d. At 
d < d"" the jamming limits decrease as the value of d in¬ 
creases. At d > d* the jamming limits increase as the 
value of d increases, where the value of d* depends on 
k. In the one-dimensional case, the results are in good 
agreement with the published analytical results [9[. The 


coverage and the jamming limits on a two-dimensional 
disordered lattice are similar to the one-dimensional case. 
The jamming limits decrease monotonically as the length 
of the line segments increases. 

Recently, the impact of defects on percolation in the 
random sequential adsorption of linear k-meis on square 
lattices was investigated for a particularly wide range of 
k-meis {k-mei lengths from 2 to 256) [T6|. Two different 
cases were studied: (1) where it was assumed that the 
initial square lattice was nonideal and some fraction, d, 
of the sites was occupied by nonconducting point defects 
(impurities); (2) where it was assumed that some frac¬ 
tion, d, of the sites in the k-meis themselves consisted of 
defects, i.e., was nonconducting, while the initial square 
lattice was perfect. 

Mixed site-bond percolation was studied for the RSA 
of k-meis on heterogeneous lattices with variable con¬ 
nectivity, 2: 0. The simulations were performed for 
k = 1-3 on a triangular lattice. Percolation phase di¬ 
agrams in terms of the percolation threshold Pc versus 
lattice connectivity, z, were obtained. For the RSA depo¬ 
sition of monomers onto a triangular lattice with defects 
preliminarily filled with /c-mers, the percolation [k = 3- 
24) [TsI and the jamming {k < 50) [lol were investigated 
by means of Monte Carlo simulations. The nonmono¬ 
tonicity of the percolation threshold as a function of the 
impurity concentration was observed [TsI . 

The RSA of polydisperse mixtures of k-meis has also 
been extensively investigated 1,0 0 . A phase di¬ 
agram separating percolating from non-percolating re¬ 
gions for mixtures of monomers and k-meis on a square 
lattice {k = 2-7) has been also determined [ 2 ^. The 
jamming coverage for a mixture has been found to be 
greater than the jamming coverage either of the compo¬ 
nents making the mixture [5[. On the other hand, the 
percolation threshold for a mixture was slightly greater 
than that of the longest k-mei 0 . The continuum ran¬ 
dom sequential adsorption of polymer on a flat and homo¬ 
geneous surface has also been studied [22| . The polymers 
were modeled as stiff or flexible chains of monomers. 

However, most of the previous studies have been de¬ 
voted to the 2D deposition of randomly oriented aniso- 
metric particles on the substrates. In several works, prob¬ 
lems with unequal probabilities for the orientation of the 
deposition of k-meis along different directions of the lat¬ 
tice have also been analyzed. The degree of anisotropy 
can be characterized by the order parameter, s, defined 
as 


where N\ and N_ are the numbers of vertically and hor¬ 
izontally oriented particles. 

The anisotropic RSA of dimers on a square lattice has 
additionally been investigated using both Monte Carlo 
simulation and time series expansion [6[. The problem 
of the anisotropic RSA of dimers on a square lattice has 
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been studied [ 2 ^. Data from the Monte Carlo simula¬ 
tions evidence that the orientational order parameter, 
5 , influences both the values of the percolation thresh¬ 
old Pc and the jamming concentration pj. In particu¬ 
lar, the value of Pj decreases as the order parameter, 
s , increases [^. An interesting finding is that in 
the limit of s ^ 1, the asymptotic fraction of dimers 
with horizontal direction does not vanish but equals 
e“^[l — exp(—2e“^)]/2 « 0.016046 [H. The properties 
of the anisotropic RSA of flexible k-meis on a 2D tri¬ 
angular lattice have been studied numerically by means 
of Monte Carlo simulations [5|. It was shown that the 
relaxation time to the jamming limit increases with the 
degree of anisotropy of the elongated particles. 

The effect of anisotropic deposition of k-meis on the 
jamming and the percolation [25| has also been in¬ 
tensively studied by means of Monte Carlo simulations. 
These detailed studies revealed that the RSA model does 
not allow preservation of the order parameter, s, and, in 
fact, the substrate ‘selects’ the k-mei with appropriate 
orientation, resulting in deviation of the predetermined 
order parameter, s, from the one actually observed one, 
sq. a special variant of relaxation for the RSA model (the 
RRSA model) with better preservation of the predeter¬ 
mined anisotropy has been developed [ 2 ^ . In the RRSA 
model, the binding of a k-mei to the adsorbing substrate 
is strong, and the k-mei has additional possibilities for 
joining the surface after an unsuccessful attempt. 

Our aim was to study the influence of anisotropy on 
the percolation threshold and the jamming concentration 
of linear k-meis onto a square lattice with point defects. 
The results for the multiple cluster relaxation random 
sequential adsorption (RRSA) model and the single clus¬ 
ter relaxation (RSC) model were compared. In the RRSA 
model, the deposition was distributed over different sites 
of the lattice, whereas in the RSC model the cluster grows 
at its perimeter, starting from an active seed in the cen¬ 
ter of the lattice. In fact, the RSC model is an extension 
of the Eden growth model [26| for the case of k-meis on 
a lattice with defects. 

The rest of the paper is organized as follows. Section HIl 
describes our RRSA and RSC models of k-mei deposi¬ 
tion. The results obtained using finite size scaling theory 
and dependencies of the percolation threshold, Pc, and 
the jamming concentration, pj, vs the order parameter, 
5 , and the defect concentration, d, are examined and dis¬ 
cussed in detail in Section Hill We summarize the results 
and conclude the paper in Section HVl 


II. METHOD 

In both the RRSA and the RSC models, the square lat¬ 
tices were initially randomly filled with point defects at a 
given concentration, d. Before deposition of the /c-mers, 
we choose an appropriate orientation in the horizontal 
or vertical directions according to the given value of the 
order parameter s. The lattice sites were then occupied 


by the addition of k-meis. The defects and previously 
deposited k-meis can inhibit the deposition of newcom¬ 
ers. If the attempt was unsuccessful, a new lattice site 
(the RRSA model) or a new empty cluster perimeter site 
(the RSC model) was randomly selected until the object 
could be deposited. The objects may move over all the 
substrate (RRSA model) or over all the cluster perime¬ 
ters (RSC model) searching for a sufficiently large empty 
space. The other specific details are presented below. 


A. RRSA model 

In the RRSA model, the deposition terminates when a 
jamming state is reached along one direction [^. We 
considered a lattice with periodic (toroidal) boundary 
conditions to eliminate the border effects and, in contrast 
to [13 , treated spiral clusters as wrapping (percolating). 
We checked the percolation in two perpendicular direc¬ 
tions and used two criteria: there is percolation in both 
directions (criterion AND), there is percolation at least 
along one direction (OR). 

Examples of the jamming states (upper row) and the 
percolation clusters (bottom row) at different values of 
the order parameter, s, are presented in Eig. [TJ The 
concentration, d, of defects on the lattice is 0.09. 

Eor each given order parameter, 5, and concentration 
of defects, d, we filled the lattice with k-meis to a con¬ 
centration p, or to jamming and checked whether there 
was percolation. We repeated this 1000 times and found 
the probability, R(p), that percolation occurs at a par¬ 
ticular concentration of the k-meis. The abscissa of the 
inflection point of the curve was treated as the estimate 
of the percolation threshold for the given lattice size. We 
used lattices of sizes L = 100/c, 200/c, 400/c and performed 
finite-size analysis to obtain the percolation threshold at 
the thermodynamic limit {L oc) Pc oc ^ where 

= 4/3 is the critical exponent of the correlation length 
for the 2d random percolation problem [T^ • All probabil¬ 
ity curves intersect each other at one point with the co¬ 
ordinates {pc^R^)^ where R* = R{pc) ~ 0.78. The value 
of R* depends on how the percolation is defined • The 
intersection point offers another way to estimate the per¬ 
colation threshold. The examples of probability curves 
and scaling are shown in Eig. [2j 

We compared our numerical simulation for completely 
aligned k-meis (5 = 1 ) with the analytical results @ 
and found that there is not visible difference between the 
analytical and numerical results. 

B. RSC model 

The Eden algorithm [26| was applied to grow a single 
cluster from an active seed /c-mer on the square lattice. 
In the RSC model, the initial seed k-mei is placed at the 
center of the lattice and it has 2k + 2 initial perimeter 
sites. The RSC algorithm uses the following steps: 




FIG. 1. Examples of the jamming states (upper row) and the percolation clusters (bottom row) of anisotropically deposited 
/c-mers on a square lattice with defects, (a) s = 0.25, (b) s = 0.5, (c) s = 1. The lattice size is 1024 x 1024, /c = 8, a fragment 
of the lattices with 128 x 128 sites is shown. The concentration of defects on the lattice is 0.09. Online: Horizontal k-meis are 
shown in red, vertical k-mers are shown in blue, k-mers not belonging to the percolation cluster are shown in the same colors 
but with a different hue, empty sites are shown in white, defects on the lattice are shown in black. Print: gray-scale. 






FIG. 2. Probability curves for /c = 4, s = 0.75. (a) R{p). (b) R{d). Insets: scaling. The statistical error is smaller than the 
marker size. 


Step 1: Randomly choose one empty cluster perimeter 
site and try to add a new k-mei to any available 
lattice sites, if such the deposition is not inhibited. 
Repeat this step using the same orientation of k- 
mer until successful attempt. 

Step 2 : Determine the new cluster perimeter sites. 

Step 3: Continue the previous steps until there remain 
no untested cluster perimeter sites or the cluster 


reaches one of the four boundaries of the lattice. 

Examples of clusters in the absence and presence of 
defects for the RSC model {k = 8) with different or¬ 
der parameters, s, are shown in Fig. [H When defects 
are absent, an infinite cluster can grow on the lattice 
and the formation of stacks of identically oriented k-meis 
(horizontal and vertical) can be observed. These stacks 
are rather similar to that observed for the RSA or the 
RRSA models [HI • Visual observations of the pat- 
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FIG. 3. Examples of the clusters in the absence (upper row) and in the presence (bottom row) of defects for the RSC model at 
different values of the order parameter: (a) s = 0.0, (b) s = 0.5, (c) s = 1. The lattice size is 1024 x 1024, /c = 8, a fragment of 
the lattices with 128 x 128 sites is shown. The concentration of defects on the lattice, d, is near the critical dc'. 0.215 (s = 0.0), 
0.186 (s = 0.5) and 0.134 (s = 1.0). Online: Horizontal k-meis are shown in red, vertical k-mers are shown in blue, empty sites 
are shown in white, defects on the lattice are shown in black. Print: gray-scale. 



FIG. 4. Probability curves Rb{d) for the RSG model, /c = 2 (a) and /c = 32 (b). The data are presented for disordered systems, 
s = 0. The statistical errors are smaller than the marker size. The lines are provided simply as eye-guides. 


terns at different concentrations of defects allow us to 
draw the conclusion that the presence of defects restricts 
the formation of stacks and can prohibit the growth of 
an infinite single cluster for some critical concentration 
of defects, dc. 

For each given anisotropy, s, and concentration of de¬ 
fects, d, we repeated the Monte Carlo experiments 1000 
times and found the probability, Rb{d), that the single 
cluster reaches one of the boundaries of the lattice at 
a given concentration of defects, d. We used lattices of 


different sizes in the interval L = 8/c-2048/c. Examples 
of the probability Rb{d) curves for a disordered system 
(s = 0) using two different values of k are shown in Fig. HI 
The values of Rb decrease as the concentration of defects 
d increases. The commonly applied procedure for de¬ 
termination of the critical concentration dc for blocking 
of the cluster growth is to use finite scaling analysis for 
a fixed value of Rb (e.g., at Rb = 0.5). However, for 
the RSC model the probability Rb{d) curve is not step¬ 
like even for infinite systems in the thermodynamic limit 
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L ^ oc. For the RSC model at L ^ oc, Rb{d) the curve 
smoothly descends as the value of d increases reflecting 
the finite probability of the blocking of cluster growth 
from the initial seed k-meis surrounded by lattice de¬ 
fects. Examples of finite-size analysis at different values 
of Rb {k = 2, s = 0) are presented in Fig. [5^1. Good 
linear dependencies of data were always observed in the 
d versus coordinates, where v = 4/3 is the criti¬ 

cal exponent of the correlation length for the 2d random 
percolation problem [a. To be definitive, we always es¬ 
timate the critical value of the concentration of defects 
dc in the limit L ^ oo at the fixed value Rb = 0.5 (see 
Fig. [IK)- At this value dc, the growth of the infinite clus¬ 
ter is suppressed with a probability of 50%. Examples of 
the dependencies for the RSC model (k = 2) at 

different values of the order parameter s {Rb = 0.5) are 
presented in Fig. |5]3. The data shows that the slope of 
the scaling can be greatly dependent on the value of s. 

The mean degree of the cluster anisotropy, d, has been 
calculated as 

S = {ry-r^)/r, (2) 

where Vy and Vx are the radii of gyration of the cluster in 
y and x directions, respectively, and r is the mean radius 
of gyration. 

The density of the cluster of /c-mers, Ps, was calculated 
inside a rectangle of size 2ry and 2rx with its center lo¬ 
cated at the geometrical center of the cluster. A finite 
scaling analysis was also carried out for the estimation of 
Ps in the limit L ^ oo. 


III. RESULTS 
A. RRSA model 

Figure [6] shows the typical dependencies of jamming 
concentration, pj^ versus the defect concentration, d, for 
different values of k and a fixed value of the order param¬ 
eter, s = 0.5. For any given k, the value of pj decreases 
as the value of d increases. The effect is more pronounced 
for the longer /c-mers, and the jamming concentration de¬ 
creases drastically with growth of k-mei length for any 
anisotropy. Similar dependencies were obtained for other 
values of the order parameter, s. 

Figure [7] shows the typical dependencies of the relative 
jamming concentration, Pj/pj{s = 0), versus the order 
parameter, s, for different defect concentrations, d, for 
k = S. For large concentrations of defects (d > 0.3), the 
jamming concentration decreases monotonically as the 
extent of anisotropy increases, whereas for small concen¬ 
trations of defects (d < 0.3), the jamming concentration 
has a minimum, near the value of s = 0.5. 

In the RRSA model, there exists a critical concentra¬ 
tion of defects, dc, above which the percolation, i.e. for¬ 
mation of an infinite cluster, is impossible even in the 
jammed state. An example of the dependencies of the 


TABLE I. Estimated parameters and a and, a in the 
power law function Eq. (3] The data are presented for RSC 
and the RRSA [U models at s = 0 and s = 1 and for the 
one-dimensional (ID) model 0 . In all cases the coefficient 
of determination, p, was greater than 0.999. 



Pj{k = oo) 

a 

a 

s 

RSC RRSA 

RSC RRSA 

RSC RRSA 

0 

0.67(6) 0.652(8) 

0.41(8) 0.417(5) 

0.77(7) 0.713(7) 

1 

0.745(0) 0.747(2) 

0.27(3) 0.235(3) 

0.93(1) 1.016(6) 

Id 

0.74759792 

0.227(7) 

1.011(1) 


percolation threshold, Pc, and jamming concentration, 
Pj, as functions of the defect concentration, d, for the 
given order parameter, s = 0.5, and k = 2 is presented 
in Fig. IHl 

For the short k-meis {k ^ 8), the percolation thresh¬ 
old, Pc, is almost insensitive to the defect concentration, 
d, for any anisotropy, s [2^ (e.g., see Fig. [8] for A: = 2). 
For long k-meis (A: > 8), Pc(d) has a maximum between 
d = 0 and d = dc [l^ - 

Figure [9] presents percolation phase diagram using as 
coordinates the critical percolation concentration, pc, 
versus the critical concentration of defects, dc, at dif¬ 
ferent values of the order parameter, s, and lengths of 
the A:-mers. The increase of order parameter, 5, is al¬ 
ways accompanied by an increase of Pc- However, the dc 
versus s behavior is different for short and long A:-mers. 
Moreover, it is interesting that for partially oriented sys¬ 
tems {s < 1) the Pc{dc) dependencies at fixed values of 
s go through a minimum at the values of k lying in the 
interval between 9 and 16. 


B. RSC model 

Preliminary investigations were performed for zero de¬ 
fect systems (d = 0). Figure [TOl presents examples of the 
degree of the cluster anisotropy, d, versus the order pa¬ 
rameter, s, for different values of k for such zero defect 
systems (d = 0). The value of 6 increases as s increases 
and the effects become more pronounced for the longer 
A:-mers. At large values of s, the shape of the clusters 
becomes ellipse-like (see inset to Fig. [TQ|) . 

Figure [TT] compares both the density of the cluster, p^, 
versus k for the RSC model and the jamming concentra¬ 
tion, pj^ versus k for the RRSA model. It is worth 
noting that the values of Ps and pj decrease as an in¬ 
verse power of the linear segment length k (see, inset to 
Fig. [TT]) and approach their limiting values p^ and p^ 
with increasing k: 

Ps,j{k) =P^j+a/k°‘. (3) 

The parameters p^ and p^, a and a are presented in 
Table H The exponent a is not universal and depends 
upon the order parameter, s. It had been suggested that 
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FIG. 5. Concentration of defects, dc, versus the scaled size of the system for the RSC model {k — 2) at different values 
of probability Rb (s = 0) (a) and at different values of the order parameter s (Rb = 0.5) (b). Here, ly = 4/3 is the critical 
exponent of the correlation length [T^. 



FIG. 6. The jamming concentration, pj, versus the defect 
concentration, d, for different values of k and a fixed order 
parameter s = 0.5. The solid lines present are simply as an 
eye-guide. The statistical error is smaller than the marker 
size. 



FIG. 7. The relative jamming concentration, pj/pj{s = 0), 
versus the order parameter, s,at different concentrations of 
defects, d. The length of the k-mevs equals 8. The solid lines 
are simply to provide eye-guidance. The statistical error is 
smaller than the marker size. 


that such power behavior indicates the presence of a frac¬ 
tal structure of the cluster networks in both the RSC and 
RRSA models with the fractal dimension df = 2 —a H. 
For completely ordered systems {s = 1), the value of df 
is close to 1 as expected for the ID problem. For disor¬ 
dered systems at s = 0, the values of a are noticeably 
smaller than 1 for both the RSC and the RRSA models 
and so the fractal dimension, d/, in this case is higher 
than 1. 

Figure [12] compares the density of the cluster, ver¬ 
sus the concentration of defects, d, for the RSC model at 
different values of the linear segment length k for disor¬ 
dered 5 = 0 (solid lines) and completely ordered 5 = 1 
(dashed lines) systems. The values of Ps decrease as the 
value of d increases and there exists some limiting con¬ 
centration of defects, dc, that suppresses the growth of 


an infinite cluster. At any given value of d, the impact of 
ordering on the value of Ps depends both upon the values 
of k and d (Fig. [12]). 

The critical concentration of defects, dc, decreases as 
the values of k and 5 increase, and approaches zero in 
the limit of very long /c-mers, k ^ oo (Fig. [T3|). It is in¬ 
teresting that a rather good linear relation between 1/dc 
and k was observed at different values of 5 : 

d~^ = d~^{k = 1) + a{k — 1), (4) 

where d-\k = 1) = 2.488 ± 0.002 and a = 0.282 ± 
0.004 (5 = 0.0), a = 0.371 ± 0.004 (5 = 0.5) and 
a = 0.794 ± 0.001 (5 = 1.0). Note, that the value 
1 — dc{k = 1) = 0.598 ± 0.003 is close to the value of the 
threshold concentration 0.5927 for ordinary monomer 
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FIG. 8. The percolation threshold, pc^ and the jamming con¬ 
centration, pj^ as a function of the defect concentration, d, for 
the anisotropy s = 0.5 and k — 2. The hatched region cor¬ 
responds to the concentrations of the objects and the defects 
when percolation is possible. Here dc is a critical concentra¬ 
tion of defects that suppresses the formation of an infinite 
(percolation) cluster. The solid lines are simply present as 
eye-guides. The statistical error is smaller than the marker 
size when not shown explicitly. 



FIG. 9. The percolation concentration, pc, versus the critical 
concentration of defects, dc, for different values of the order 
parameter, s, and lengths of the /c-mers. The lines are simply 
present as eye-guides. The statistical error is smaller than the 
marker size. 


percolation [a. 

The density of the clusters, ps, at the critical concen¬ 
tration of defects, dc, for the RSC model depends in a 
complex manner on s and k (Fig. [14]). For disordered 
systems (s = 0), the value of Ps tends to zero in the lim¬ 
its of the very long /c-mers, k ^ oo and very small critical 
concentrations dc ^ 0. The formation of such ‘empty’ 
or loose clusters can be explained by analyzing the clus¬ 
ter patterns presented in Fig. |3l In zero defect systems 
(d = 0), compact clusters form with significant stacks 
of identically oriented k-meis (horizontal and vertical). 



0.0 0.2 0.4 0.6 0.8 1.0 

S 


FIG. 10. Examples of the degree of the cluster anisotropy, 
d, versus the order parameter, s, for different values of k for 
zero defect lattices (d = 0). Inset: examples of the clusters 
for /c = 8. 



FIG. 11. Density of cluster, ps, versus the value of k for the 
RSG model, and the jamming concentration, pj, versus the 
value of k for the RRSA model [Ij. The data are presented 
for disordered s = 0 and completely ordered s = 1 zero defect 
systems (d = 0). Inset: ps — pT Pj ~ pT v^i'sus k. Here, 
pj° and pY are the limiting values at k ^ oo. 

These stacks are responsible for the high values of p^. 
The introduction of defects strongly prevents the forma¬ 
tion of stacks (see, e.g., Fig.|3^) and loose networks with 
near zero density are formed. On the other hand, denser 
clusters are formed for ordered systems with ps ~ 0.065 
at s = 0.5 and Ps ~ 0.38 at s = 1.0. 
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FIG. 12. Density of cluster, pc, versus the concentration of 
defects, d, for the RSC model at different values of the lin¬ 
ear segment length, k, for disordered s = 0 (solid lines) and 
completely ordered s = 1 (dashed lines) systems. 



FIG. 14. Density of clusters, ps, versus the critical concen¬ 
tration of defects, dc, for the RSG model at different values 
of the order parameter, s, and of k. Arrows show the values 
of ps in the limit of very long /c-mers, k ^ oo. 



FIG. 13. Gritical concentration of defects dc versus the value 
of k for the RSG model at different values of the order pa¬ 
rameter, s. Inset: 1/dc versus k dependencies. 
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IV. FINAL REMARKS AND CONCLUSIONS 

For the problem of k-meis deposition, the RRSA and 
the RSC models have huge differences in their deposition 
rules and structure of the clusters they generate. In the 
RRSA model, random deposition of k-meis is performed, 
with multiple clusters being formed that can consolidate 
in the course of their growth. Intensive studies with 
the RRSA have shown that the jamming concentration 
continuously decreases as the lengths of the k-meis in¬ 
crease [ 2 ^ . The percolation threshold initially decreases 
and then increases with increasing values of k [25[ . For a 
completely disordered system, i.e. at 5 = 0, a conjecture 
has been offered that the formation of an infinite cluster 
is impossible when k exceeds approximately 1.2 x 10^ [25[. 
On the other hand, in the RSC model the single cluster 
grow at its perimeter by an Eden growth process [26[ and 
the formation of an infinite cluster is not restricted for 
any values of k. This model assumes strong attractions 
of the newcomer to the previously deposited particles at 
their lateral boundaries. The Eden-like clusters display 
noticeable shape anisotropy at high values of the order 
parameter, s. 

It is amazing that, in spite of these differences between 
the RRSA and the RSC models, they display some gen¬ 
eral similarities. In the absence of defects on the sub¬ 
strate (for zero defect systems, d = 0), the jamming 
concentration, pj, for the RRSA model and density of 
the cluster, ps, for the RSC model, display similar de¬ 
pendencies on the linear segment length, /c, and on the 
order parameters (see, Eig. [I2|. The impacts of defects 


on the percolation and jamming characteristics are also 
fairly similar. Eor both models the suppression of the 
growth of an infinite (percolation) cluster at some criti¬ 
cal concentration of defects, dc, can be observed. Phase 
diagrams in the form of the Pc(dc) (RRSA model) and 
Ps(dc) dependencies for different values of s and k were 
determined. 

Eor the RRSA model, the value of Pc ranges in the in¬ 
terval 0.465 — 0.58 {k = 2 — 32) and the value of dc 
decreases for short k-meis {k < 16) but increases for long 
k-meis {k > 16) as the value of s increases. Moreover, for 
short /c-mers, the percolation threshold is almost insen¬ 
sitive to the defect concentration for any values of s. Eor 
the RSC model the value of ps at the critical concentra¬ 
tion of defects, dc, depends in a complex manner on the 
values of s and k. Eor disordered systems (5 = 0 ), the 
value of Ps tends to zero in the limits of very long /c-mers, 
k ^ 00 . This reflects a suppression of k-mei stacking by 
the defects that results in the formation of loose clusters 
with very low density. 
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